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Abstract 

In this paper we introduce a new model of data packet 
transport, based on a stochastic approach with the 
aim of characterizing the load distribution on com- 
plex networks. Moreover we analyze the load stan- 
dard deviation as an index of uniformity of the dis- 
tribution of packets within the network, to character- 
ize the effects of the network topology. We measure 
such index on the model proposed by Watts and Stro- 
gatz as the redirection probability is increased. We 
find that the uniformity of the load spread is maxi- 
mized in the intermediate region, at which the small 
world effect is observed and both global and local ef- 
ficiency are high. Moreover we analyze the relation- 
ship between load centrality and degree centrality as 
an approximate measure of the load at the edges. 
Analogous results are obtained for the load variance 
computed at the edges as well as at the vertices. 



1 Introduction 

Congestion in real networks is a complex phe- 
nomenon that depends on a large number of vari- 
ables. Here we are interested in understanding how 
the underlying structure of the network can itself 
have an influence on congestion. The network is sup- 
posed to be crossed by information units or packets. 
We will consider a packet transport model on a given 
graph and will study the effects of changes in the 
graph topology on the packets distribution. With 
the aim of understanding the effects of the network 
topology on congestion, we will keep constant other 
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variables as the rate at which packets are generated, 
the distribution of sources and destinations within 
the network, and the routing algorithm to go from 
the former to the latter. We will show that each ver- 
tex can draw toward itself an higher or lower flow 
of packets, according to its own position within the 
graph. 

2 The network as a Markov 
Chain 

In the recent literature on networks several dynamic 
packet transport models were analyzed with the aim 
of describing realistic communication over the net- 
work Q21 00 EU On the other hand in © 
some static parameters were introduced in order to 
characterise the load distribution over the network 
without considering explicitly its dynamics. 

We will suggest here an alternative and interme- 
diate packet transport model based on a stochastic 
approach. 

We consider a network with N vertices. We make 
the hypothesis that the origin and the destination 
of each packet are chosen with uniform probability 
within all the vertices in the network. We suppose 
that at each time step, each vertex generates N — 
1 packets, each one addressed to every other node 
in the network so that the total number of packets 
generated at each time step is equal to N(N — 1). 

Then packets are supposed to be routed from origin 
to destination along the geodesic, i.e. according to a 
minimum distance criterion. We suppose that pack- 
ets are delivered according to a dynamic process that 
takes into account the actual time needed for packets 
to travel across the network. We consider the case 
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of a parallel process, in which at each time unit all 
packets move simultaneously, every one crossing one 
edge. We assume that the time in which a packet 
goes through exactly one edge is equal to one time 
unit. Specifically we assume that at each time step: 

1. Packets are generated at every vertex. 

2. Packets in a given node but not delivered yet to 
their final destinations, are routed to the neigh- 
bors of that node that are nearest to the desti- 
nation. 

3. Packets delivered to their final destination are 
removed from the network. 

In 0, the load at a given vertex v, £(v), was de- 
fined as the number of shortest paths between pairs 
of nodes crossing it. Analogously we will define here 
the load at the edges £(e), as the number of shortest 
paths crossing a given edge. 

Let J(i) be the subgraph of the vertices j adja- 
cent to i. As a packet is crossing i, it can take one 
of the directed edges ey, j E J(i). We assume that 
each edge in the graph is crossed with frequency 
proportional to its load /(e). Consistently with that 
assumption, we will estimate the conditional proba- 
bility of a packet starting from i, going through one 
of the edges ey, as: 



Using QJ, we can now associate to the network of 
interest a Markov Chain; each node being represented 
as a state of the chain with the probability of going 
from a state i to a state j being given by p(eij/i) for 

jeJ(i). 

The Markov Chain represents a dynamic model de- 
scribing the evolution of the packets distribution in 
time. We can obtain numerically the limit distribu- 
tion A, i.e. the probability of being asymptotically 
in a certain state/node of the Markov Chain. We 
wish to emphasize that A can be used as an alterna- 
tive measure of the load at vertices. In particular, we 
found A to differ from the previous estimation of the 
load at vertices [S]. 



We can give a physical explanation of that, accord- 
ing to the parallel nature of the packets exchange over 
the network. 




Figure 1 : A parallel model of packets exchange in a sam- 
ple network of three nodes. We represent with the same 
color packets generated at a given vertex. In the left panel 
we represent the distribution of packets after their gen- 
eration. It is worth noticing that the distribution does 
not change if we look at it after packets have been deliv- 
ered. As we can see, packets just having been generated 
and the ones delivered to destination are never present 
simultaneously in the same node at any given time. 



Fig.l shows the load distribution of packets in 
a representative network on which packets are ex- 
changed in a parallel way. We observe that for each 
packet travelling from a vertex j to a vertex j, there 
is another one travelling from j to i. This implies 
that packets just having been generated and the ones 
delivered to destination are never present simultane- 
ously in the same node at any given time. 

Let's define A„ the total load at the vertices in the 
network, computed as the sum of the loads over all 
the vertices within the network, i.e. A v = ^2 v £(y). 
We have that: 

A v = d- N(N - 1), (2) 

where d is the characteristic path length associated to 
the graph. From (J2J , we have that the average load at 
vertices is given by I = d-(N-l). We observe that the 
average load is proportional to the characteristic path 
length of the graph: shortening mean distances be- 
tween vertices implies that packets travel from origin 
to destination in fewer steps. Hence, a load reduction 
is observed on the whole network. 
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Moreover, it is worth mentioning here that expres- 
sion (0) is identical to Little's Law ^U] m queueing 
theory, i.e. 

V = t-N(N -1), (3) 

where the total number of packets in the network V 
is replaced with A v and r, the mean time spent by 
a packet going from origin to destination in Little's 
law, is replaced with d, under the hypothesis that 
each edge is crossed within the same amount of time. 

Note that, according to the newly introduced load 
definition, we find that, contrary to what presented 
in [5], the total load at the vertices is equal to the 
total load at the edges A e = J2 e l(e), i.e. A e = A„. 
This is expected as in our model, the total amount 
of packets within the network at any time has to be 
the same. 

3 Uniformity of load distribu- 
tion in small world networks 

We will try now to assess whether the network topol- 
ogy can itself cause a more or less uniform distri- 
bution of packets in the network. In [5] the ef- 
fects of variations of load distribution in random and 
scale free networks were analyzed. The former are 
characterized by a Poisson degree distribution while 
the second have a degree distribution of the form 
P(k) ~ /c~ 7 . From a communication point of view, it 
would be desirable to have a uniform load distribu- 
tion in order to exploit evenly the network resources 
(nodes and edges). In particular we will try to un- 
derstand how uniform over a given network the load 
distribution is, because of its topology. 

Specifically, an high variance of that distribution 
should indicate an unfair use of the network, and 
could therefore point out a possible cause of conges- 
tion. 

In PP, the maximum load £ ma x was proposed as a 
main index for characterizing the network structure 
and an high value of this parameter was claimed re- 
sponsible for congestion. Nevertheless, this parame- 
ter does not describe the whole structure of the graph, 
giving information only on the load level at the hub. 
For that reason here we will evaluate the load stan- 



dard deviation as defined in that refers to the 
network as a whole. In order to make it insensitive 
to the average values of i, we will evaluate cr{€), the 
standard deviation of the load appropriately normal- 
ized with respect to its mean. 

The standard deviation of the load distribution 
<j(£) was shown to have an effect on network per- 
formances in terms of throughput and delivery time. 
Namely an high value of the load variance, that is 
typical of scale-free networks was seen to be able to 
worsen strongly the network performances. 

We analyze here the variation of load distribution 
both at the edges and at the vertices due to changes 
in the underlying topology. 

A common feature of most real networks is that 
they are small world, i.e. the average distance be- 
tween randomly chosen nodes is generally low and 
increases only logarithmically with the network size. 
Moreover real networks show high clustering or tran- 
sitivity, that is, high probability that two randomly 
selected neighbors of a node are also neighbors of each 
other. Watts and Strogatz (WS) 0] showed that, 
starting from a regular nearest neighbor circle, a few 
random redirection were sufficient to get the small 
world effect, without compromising the clustering. 

We will repeat the experiment in .4], measuring 
also the normalized standard deviation of the vertices 
load, a(£(v)) and the edges load a(£(e)). 

In Fig. 2 the average load in the graph is shown 
as the redirection probability p is varied. It can be 
observed that the average load behaves as the char- 
acteristic path length of the graph, decreasing as p 
increases. 

The normalized standard deviation of the vertices 
load is shown in Fig. 3. Here, we see that such quan- 
tity increases at first, then over a certain threshold 
value of the redirection probability, it is progressively 
reduced. In fact, as the redirected edges are few, the 
global reduction of the characteristic path length is 
achieved through the exploitation of the same short- 
cuts and the vertices situated at their endpoints. 

A further increase of p (i.e. of the number of redi- 
rected edges) causes an increase of the number of 
different paths linking pairs of nodes in fewer steps. 
Thus, for higher p, an higher load uniformity (lower 
a(£(v)) is observed. In other words we observe that, 
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Figure 2: The WS model. A regular lattice with 1000 
vertices and 10 links per vertex is rewired with probability 
p. Average vertices load £ = d ■ (N — 1) is reported as a 
function of p. The results are obtained over 100 diverse 
experiments. 
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Figure 3: Normalized standard deviation of vertices load 
a(£(v)) is reported as varying the redirection probability, 
under the same hypothesis of figure 2 



starting from a regular configuration, and adding pro- 
gressively some disorder to it, the parameter increases 
at first and then decreases after a critical point. So 
we have that at intermediate p, at which local and 
global efficiency are both maximized we S e t the 
least uniform load distribution between vertices. In 
Fig. 4 we report the median (in blue) of the load dis- 
tribution with the 25th and 75th percentile (in red), 
showing that the interquartile range of the sample is 
larger for intermediate values of p. 

We can argue thus that the small-world effect, i.e. 
the reduction of the average distance between vertices 
caused by the addition of a relatively small number of 
shortcuts, can be associated with an unfair spread of 
the load distribution across the network. This would 
cause an unfair exploitation of some of the network 
resources, making some areas of the network more 
likely to be interested by higher traffic and therefore 
by congestion phenomena. 

As shown in Fig. 5, similar results are obtained 
when the load standard deviation at the edges (rather 
than at the vertices) is computed. 

4 Load Distribution And De- 
gree 

As discussed above, by adding shortcuts at a given 
vertex, a load increase is observed at the node as 
well as at the edges leaving from it. So we can notice 
two different effects of rewiring links on the load at 
the edges: (i) a bigger amount of packets is routed 
by the starting node through the outgoing edges; (ii) 
the load is shared between an higher number of edges. 
Therefore, the load distribution at the edges is influ- 
enced by both the load and the degree at the vertices 
which they leave from. 

From this point of view, it is worth understanding 
the relationship between the two main measures of 
point centrality j the one based on the degree and 
the other on the betweenness. 

A first contribution in that sense was given in [8], 
where it was shown that, in the case of a power law 
degree distribution, the load scaling were character- 
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Figure 4: The median of load distribution is reported 
as varying the redirection probability (blue line), under 
the same hypothesis of figure 2. Upper and lower red 
circles represent the 25th and the 75th percentile, with 
the length of the red segments indicating the interquartile 
range. 
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Figure 5: Normalized STD of edges load a(£(e)) is re- 
ported as varying the redirection probability, under the 
same hypothesis of figure 2 



ized by a function of the degree k: 

P(£) oc k v , (4) 

This means that nodes with more incident edges are 
the ones that draw on themselves more load, accord- 
ing to an assigned power of the degree. 

Here, to define a more general link between load 
and degree, we introduce a new index, the load-degree 
ratio r(i) at each vertex i, defined as 

r[l) ~W)~ W) ' ( ' 

This parameter measures the number of packets at 
each vertex, with respect to the number of incident 
links and moreover it is an average value of the load 
at the edges incoming or outgoing from i. 

We observe that, when the redirection procedure 
is applied, the degree of some vertices increases (as 
new edges are connected to them). Correspondingly 
a significant increase of the load at vertices is also 
detected. Therefore, the load-degree ratio of those 
vertices increases as well. On average this means that 
the load at the edges leaving from overloaded vertices 
becomes larger. This explains the behavior of a(£(e)) 
in Fig. 5 below the threshold value of p. 

On the contrary, for further variations of the redi- 
rection probability, as the network approaches a 
random configuration, a(£(e)) decreases similarly to 
what observed for a(£(v)). 

5 CONCLUSION 

In this work we have introduced a new load definition 
more suitable in the case of a dynamical process in 
which packets travel across the network in parallel. 
This new formulation differs from previous ones as it 
takes into account in the computation alternatively 
the contribution of the packets yet generated or of 
the ones delivered to destination; in fact in a parallel 
process packets outgoing from a given vertex move at 
the same time as the incoming ones arrive. 

Moreover we have introduced two new indices with 
the aim of evaluating the role of the network topol- 
ogy in influencing the packets distribution within the 
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network, the standard deviation of the load at the 
vertices a(£(v)), and at the edges a(£(e)). We have 
reworked through the WS experiment, computing the 
new quantities. We have noticed that as the small- 
world effect can be obtained by the addition of a rel- 
atively low number of shortcuts, this is made at the 
expenses of the uniformity of the load distribution 
within the network. 
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